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The domain wall fermion formulation exhibits full chiral symmetry for finite lattice spacing except for the effects 
of mixing between the domain walls. Close to the continuum limit these symmetry breaking effects should be 
described by a single residual mass. We determine this mass from the conservation law obeyed by the conserved 
axial current in quenched simulations with j3 = 5.7 and 6.0 and domain wall separations varying between 12 and 
48 on 8 3 x 32 and 16 3 x 32 lattices. Using the resulting values for the residual mass we perform two complete and 
independent calculations of the pion decay constant. Good agreement is found between these two methods and 
with experiment. 



1. INTRODUCTION 

An important feature of both domain wall and 
overlap fermion formulations is that the chiral 
limit is disentangled from the continuum limit. 
For the case of the domain wall fermions, the chi- 
ral symmetry breaking is solely controlled by the 
amount of coupling between the domain walls [jlj. 
To exploit the chiral properties of the domain wall 
fermions, it is essential to understand these chiral 
symmetry breaking effects. 

We reported earlier (||| that for quenched sim- 
ulations, although the quantity m^.(m/ = 0) ob- 
tained from the pseudoscalar density correlator 
decreased as the separation between the domain 
walls L s was increased, it did not appear to van- 
ish in the L s — > oo limit. Further investigation 
shows that the chiral limit of the pion mass is 
distorted by many factors. Besides the order 
a 2 effects and finite volume effects common to 
both dynamical and quenched simulations using 
domain wall fermions, the presence of topologi- 
cal near-zero modes in the quenched approxima- 
tion significantly affects the determination of the 
pion mass m] . These zero-mode effects can be 
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suppressed by going to a larger volume. Other 
non-linear behaviors caused by the absence of 
the fermion determinant in quenched simulations, 
such as "quenched chiral logs" , may also play an 
important role in the behavior of the pion mass 
as the chiral limit is approached. These factors 
make it difficult to determine the degree of chi- 
ral symmetry breaking from examination of the 
pion mass in the quenched chiral limit. The tech- 
nique described in detail below quantifies the chi- 
ral symmetry breaking effects more precisely by 
a single residual mass measured from the extra 
term in the divergence of the axial current || . 

2. RESIDUAL CHIRAL SYMMETRY 
BREAKING 

Close to the continuum, the effect of mixing 
between the domain walls produces chiral sym- 
metry breaking terms with a coefficient m res , the 
residual mass, in the low energy effective La- 
grangian. This results in an effective quark mass 
of m c ff = nif + ?Tircs entering all low momentum 
Green's functions. 

We can measure m rcs starting from the diver- 
gence of the axial current, given by|6| 

A^(x) = 2m f J«(x)+2J« q (x). (1) 

Compared with the corresponding continuum ex- 
pression, the additional term from Jg q , referred as 
the "mid-point" contribution, should be equiva- 
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lent to m les J^ up to 0(a 2 ). Therefore when t is 
large enough such that low energy physics domi- 
nates, the ratio 



R(t) 



(E s ^(f,^°(o)) 



(2) 



should be equal to m rcs . Fig. [I] shows that the 
plateau for R(t) at (3 = 6.0, L s = 16 starts from 
t = 2, and at (3 — 5.7, L s ~ 48 it starts from 
t = 4. 
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Figure 1. versus t at /3 = 6.0, L s = 16 

(circles) and (3 = 5.7, L s = 48 (diamonds). 

Being part of the mass term of the low en- 
ergy effective Lagrangian close to the continuum, 
m res should provide a universal description of the 
residual chiral symmetry breaking effects in all 
long distance observables. This is demonstrated 
in Figure || where m YOS shows little mj or volume 
dependence at (3 = 5.70, L s = 48. 

Figure || shows the L s dependence of m rcs at 
(3 = 6.0. The residual mass decreases with L s , 
but the data points are poorly fit by a simple 
exponential. If we add a constant to the function, 
a much better fit can be obtained. This constant, 
which is the residual mass at infinite L s , is about 
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Figure 2. The residual mass m res versus m/ at 
f3 = 5.70, L s — 48. The open circles are obtained 
from 106 configurations with volume 16 3 x 32, 
while the closed diamonds are obtained from 335 
configurations with volume 8 3 x 32. 



1 MeV. The five points can also be fit extremely 
well to a double-exponential function; this may 
suggest the existence of multiple decay modes. 
The exact asymptotic form for function 
of L s still needs further investigation. However 
for L s = 16, m ICS is already as low as 3.87(16) 
MeV in the MS scheme at 2 GeV, about 1/30 of 
the strange quark mass. This value is so small 
that the residual chiral symmetry breaking is not 
very important for these simulations. Even at the 
stronger coupling (3 = 5.7 (a -1 ~ 1 GeV), when 
L s is increased from 32 to 48, 77i ros drops from 
0.0105(2) to 0.00688(13), which is about 1/14 of 
the strange quark mass. Good chiral properties 
can be obtained by simulating at this value of L s . 
Studies of this issue can also be found in . 

3. CALCULATION OF U 

As a test of these measurements and chiral 
properties, we calculate the pion decay constant 
using the pseudoscalar density correlator, which 
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Figure 3. The residual mass m Tes versus L s at 
j3 = 6.0. The dotted line is a simple exponen- 
tial fit. The long dashed line is an exponential 
plus constant fit. The solid line is a double- 
exponential fit. 

requires knowledge of the residual mass, and com- 
pare the results with those obtained from the ax- 
ial vector current correlator. 

Following the definition of the pion decay con- 
stant, fa can be related to the amplitude of the 
axial current correlator Aaa- 



Aaa 



i fir \2 m K 

[ Za~' 2 ' 



(3) 



In Eq. ||, the renormalization factor Za is needed 
to relate the local axial current to the conserved 
current in the domain wall fermion formulation. 
Za can be measured by studying the large t be- 
havior of the quantity Za{^), defined as the ratio 
of the coupling between the conserved current and 
the pion to the coupling between the local current 
and the pion. Figure |] shows that nice plateaus 
for Zjy(t) can be found over the ranges 4 < t < 14 
and 18 < t < 28 for both (3 = 6.0, L s = 16 with 
16 3 x 32 volume and (3 = 5.7, L s = 48 with 8 3 x 32 
volume. The L s dependence of Za at (3 — 6.0 is 
presented in Fig. || . For L s ranging from 12 



to 48, the change in Za is less than 1%. There- 
fore, the effect of finite L s on the result for Za is 
negligible. 
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Figure 4. Za{£) versus t at j3 = 6.0, L s — 16 on 
16 3 x 32 lattices (circles) and (3 = 5.7, L s = 48 
on 8 3 x 32 lattices (diamonds). 

Another method to calculate f v uses the rela- 
tion between / w and the amplitude of the pseu- 
doscalar correlator App: 



App = -(- 



-) 



(4) 



to/ + m r 

which directly depends upon the residual mass 

TO rcs . 

Results for j3 = 5.7, L s = 48 with lattice 
volume 8 3 x 32 show a discrepancy in the rrif 
dependence of f„ between these two methods. 
While the difference appears to become somewhat 
smaller at the larger volume of 16 3 x 32, we at- 
tribute the discrepancy between these two meth- 
ods of calculating f„ to 0(a 2 ) errors. 

Figure [| shows the comparison of these two 
independent calculations at the weaker coupling 
/3 = 6.0. Consistent results are seen at each va- 
lence quark mass, and the values measured by 
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Figure 5. Za versus L s at (3 — 6.0 on 16 3 x 32 
lattices. The dashed line fits the points to the 
form of A + Bcxp(— aL s ). 



extrapolation to m/ = — m rcs agree well with the 
physical value of 130 MeV, where m p is used to 
set the scale. It is important to note that these 
good properties disappear if the value for m res in 
Eq. |J is replaced by or by the x intercept of the 
simple linear fit to . Therefore, the good agree- 
ment seen between the two approaches serves as 
a consistency check of the residual mass analysis 
and a demonstration of the good chiral properties 
of the domain wall fermion formulation. 

4. CONCLUSIONS 

The chiral symmetry breaking effects in the do- 
main wall fermion formulation can be character- 
ized by a residual mass, which enters the effec- 
tive quark mass in the low energy effective La- 
grangian. This quantity can be determined accu- 
rately from the additional term in the divergence 
of axial current. Good chiral properties can be 
achieved for the quenched domain wall formula- 
tion for lattice spacing a -1 ~ 1 — 2 GeV. The 
results for further checked by our deter- 

minations of the pion decay constant. 



Figure 6. f ni K versus to/ at (3 = 6.0 on 16 3 x 
32 lattices from (Aq(x)Aq(O)) correlator (o) and 
(7r a (x)7r a (0)) correlator (o). The solid symbols 
represent the extrapolations to to/ — — m rcs and 
interpolations to the kaon mass. The diamonds 
are slightly shifted to the right for clarity. 
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